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Abstract 

Range expansion of spreading organisms has been found to follow three types: (i) linear expansion with a constant rate of 
spread; (ii) bi-phase expansion with a faster linear expansion following a slower linear expansion; and (iii) accelerating 
expansion with a continuously increasing rate of spread. To date, no overarching formula exists that can be applied to all 
three types of range expansion. We investigated how propagule pressure, i.e., the initial number of individuals and their 
composition in terms of dispersal ability, affects the spread of a population. A system of integrodifference equations was 
then used to model the spatiotemporal dynamics of the population. We studied the dynamics of dispersal ability as well as 
the instantaneous and asymptotic rate of spread. We found that individuals with different dispersal abilities were spatially 
sorted with the stronger dispersers situated at the expanding range front, causing the velocity of expansion to accelerate. 
The instantaneous rate of spread was found to be fully determined by the growth and dispersal abilities of the population at 
the advancing edge of the invasion. We derived a formula for the asymptotic rate of spread under different scenarios of 
propagule pressure. The results suggest that data collected from the core of the invasion may underestimate the spreading 
rate of the population. Aside from better managing of invasive species, the derived formula could conceivably also be 
applied to conservation management of relocated, endangered or extra-limital species. 
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Introduction 

The ability to estimate the rate of spread of an invasive species is 
important for tlie success of its management and control [1]. Early 
theories suggested that the velocity at which a species expands its 
range depends on the population growth and dispersal rates [2]. 
Models based on partial differential equations, specifically the 
reaction-diffusion (RD) model, assume a normal distribution of 
species' dispersal distances (i.e. dispersal kernel) and yield a widely- 
used formula which depicts a constant rate of spread {c = l\/rD, 
where r and D denote the intrinsic growth and diffusion rates, 
respectively) [2-4] . However, a growing body of evidence suggests 
that the rate of spread for most species may not be constant [5] . 
Shigesada et al. [6] group patterns of range expansion into three 
types: type I, linear expansion with a constant rate of spread; type 
II, bi-phase expansion with a faster linear expansion following a 
slower linear expansion; type III, accelerating expansion with a 
continuously increasing rate of spread. 

To fully comprehend the accelerating nature of type II and III 
range expansion, different dispersal strategies have been incorpo- 
rated into RD models. In particular, fat-taUed dispersal kernels (i.e. 
movements with a substantial portion of long-distance dispersal) 
have been shown to be capable of boosting the range expansion 
and are, thus, an appropriate mechanism for explaining the 
accelerating range expansion [7]. However, this explanation 



suffers from two pitfalls [8]. First, the rate of spread predicted 
from a fat-tailed dispersal kernel will keep increasing without an 
upper bound, an obvious exaggeration of the reality. Second, 
estimating the parameters of a fat-tailed dispersal kernel is difficult 
due to the obvious rarity of long-distance dispersal events [9] and 
often requires a substantial amount of recapturing records [10]. 
Clark et al. [11] tackle the conundrum of type II biphasic and type 
III accelerating range expansion by using a combined dispersal 
kernel, with the individual having a probability of p to move a 
short distance and a probability of 1 —p to move a long distance 
(see also [6,12]). This combined dispersal kernel can lead to a 
budding pattern of stratified range expansion, with the expansion 
speeding up when the buds of founding populations merge into a 
super colony. Such stratified dispersal is especially common in 
species with multiple stages (e.g. species with dispersal polymor- 
phism or different dispersal strategies at different instars; [13]) or 
with multiple dispersal vectors (e.g. combined wind and animal- 
born seed dispersal; [14]). Introduced species often experience 
stratified dispersal due to the additional human-facilitated 
translocation [15-19]. 

Two recent insights from invasion biology sketch a new 
alternative concept that could explain the range expansion that 
accelerates to a limited speed. First, propagule size (i.e., the 
number of individuals released into an introduced area) has been 
identified as one key factor of invasion success [20-23]. A large 
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propagule size can efficiently counteract the positive density 
dependence caused by the Allee and founder effects that hamper 
the establishment of initial propagules in a novel ecosystem [24— 
26]. More importandy, studies show that the initial propagule 
often consists of a suite of individuals with different performance 
ability [25,27-29], and assuming propagules with identical traits 
often leads to an underestimation of the spreading rate in animals 
[30]. It is, thus, more reasonable to conceptualize the initial 
propagule as a group of individuals with differences in their life- 
history traits. 

Second, dispersal strategy is a density- and context-dependend 
adaptive trait. Overwhelming cases support a positive density- 
dependent dispersal [31], and we expect to see a higher dispersal 
rate at the core high-density population than at the low-density 
marginal population. The interplay of local adaptation and 
environmental gradient can further lead to a context-dependent 
dispersal [32-34], with dispersal strategy and spreading behaviour 
highly sensitive to the spatiotemporal variability of habitat quality, 
especially during range shifts [35-38]. Besides the commonness of 
the density- and context-dependent dispersal strategy, in many 
invasive species, the advancing range front poses an additional 
selection force of spatial gradient, only existing temporarily while 
the range expansion is ongoing. Individuals will be sorted along 
this non-equilibrium spatial gradient according to their dispersal 
abilities, with individuals having stronger dispersal abilities more 
likely to locate at the advancing range front [39,40]. Examples of 
spatial sorting at the advancing range front are accumulating in 
literature, such as developing longer legs in cane toads (Bufo 
marinus) in Australia [41]; changing wing shape of Indian mynas 
{Acridotheres tristis) in South Africa [42,43], and more long- 
winged morphs of bush crickets in UK [44]. The range expansion 
could be accelerating due to a dynamic dispersal kernel driven by 
this process of spatial sorting for stronger dispersers at the range 
front. Indeed, based on recapturing records and metrics of spatial 
genetics, Hui et al. [10] and Berthouly-Salazar et al. [45] further 
demonstrate a changing dispersal kernel through spatial sorting in 
the invasive European starling [Stumus vulgaris) to support its 
accelerating range expansion in southern Africa. 

Given these two insights from invasion biology, the RD model 
needs to be revisited and the formula c = 2\/rD revised so that the 
velocity of accelerating range expansion can be accurately 
predicted. Here, we present a mathematical model that uses 
integrodifiFerence equations to incorporate individuals with differ- 
ent dispersal abilities in the initial propagule. The new formulae of 
the instantaneous and asymptotic rates of spread derived from this 
model include not only rates of growth and dispersal as in the 
formula for lin(-ar expansion, but also parameters depicting the 
propagule size, its composition and the process of spatial sorting. 
The classic formula c = IsfrD is shown to be a special case of the 
new ones. We further develop numerical simulations to test the 
performance of these formulae and advocate the use of them in the 
modelling and risk assessment of invasive species [46], as well as in 
forecasting the range shift of species in response to environmental 
changi;s [47,48]. 

The Model 

For simplicity, we consider the invasion of a one-dimensional 
habitat by the mixed propagule which consists of n types of 
individuals with different dispersal abilities. Let Mi(x,i) (? = 1,..., m) 
denote the population size of type i individuals at location x and 
time I. The dispersal of type i individuals is depicted by the 
dispersal kernel ki{x,y) (i.e. the probability that a type i individual 
moves from location 3) to x during a time step [7]). Dispersal can be 



density dependent [49] and sensitive to habitat quality [10], 
complicating the formulation of the spreading rate. Consequently, 
we assume homogenous habitat and density-independent dispers- 
al; that is, the dispersal kernel A, depends only on the distance 
between the locations (z= |x— Let df represent the variance of 
the dispersal kernel A,(z). For a Gaussian kernel, we have 

ki{z) = exp( — Q.Sz^/df)/^! Indf , and for a Laplace kernel, 
lid}. Hereafter d} will be referred to 



fc,(z) = exp(-zy2/^)/^ 

as the dispersal ability of type i individuals. 

The recruitment of type i individuals is governed by a non- 
negative function of population change rate gi{u\,...,u„). We 
assume that the dispersal ability is inheritable as in an asexual 
population [50,51]; that is g',(Mi,..., «„) = 0 if m, = 0. We further 
assume that the population does not suffer from the AUee effect 
and that the population change rate thus reflect the negative 
density dependence; that is, g,(ii\,..., u„)^RjUi, where Ri 
( = 8gi(0,...,0)/dui) is the population growth rate. We also assume 
that there is no trade-off l)C'twecn the dispersal ability and the 
population growth rate; that is, individuals with different dispersal 
abilities have an equal population growth rate, Ri = ...=Rn=R. 
Examples of such population change functions include the Richer 
model [52], g,(Mi,...,t(„) = M,e''^"^^, where u is the total population 
(2m,), and r ( = InR) the intrinsic population growth rate. 

Given the population size Ui{x,t) at time t, the population size at 
time t+X can be depicted by the following integrodifference 
equation [7], 

Ui{x,t+\)= ki{x-y)gi(u]_,...,u„)dy,i=\,...,n. (1) 
The total population is governed by 



u(x,t+ 1) = 



Kx,y; t)g(u)dy. 



(2) 



where k{x,y;t) {='^ki{x—y)pi{y,t)) is the expected dispersal 
1=1 

kernel at location y and time t, and Pi(x,t) the proportion type i 
individuals in the mixed population at location x and time t, 
Ui(x,t)/u(x,t). For the Ricker model, we have g{u) = ue''^"l^ . 
Evidentiy, the expected dispersal kernel not only depends on the 
dispersed distance (|x-3)|) but also the population composition 
pi(y,t) where the individuals are dispersed from, suggesting a 
dynamic dispersal kernel [8]. 

Let the initial propagule consist of Uq individuals. We consider 
two specific scenarios of mixing individuals with different dispersal 
abilities in the initial propagules. The first scenario assumes that 
the initial propagule only consists of two types of individuals, with 
their dispersal abilities d\ and d^ (i/j <rf|). More specifically, we 
assume that the majority of the initial propagule are type 1 
individuals with the dispersal ability d^, and only a small 
proportion of individuals have the dispersal ability i/fi that is, 
Mi(0,0)-I-M2(0,0)= [/o, M2(0,0)< <Mi(0,0) and 

Mi(x,0) = K2(x,0) = 0 for x^O. 

The second scenario assumes individuals with many more 
dispersal al)ilitics in the initial propagule. As the log-normal 
distribution has been confirmed for many species- and trait-level 
frequency distributions in ecological communities [53,54], we 
generated the dispersal abilities of the initial propagules as follows. 
First, for each individual, a dispersal ability d^ was randomly 
drawn from a log-normal distribution h\N{^,a) (hereafter, the 
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Figure 1. Spread of a population with two dispersal abilities. Parameter values are r = 0.15,clj ^\,d^ ^5. A-C: Population size at different 
time with the initial populations u\ (0,0) = 0.99,W2(0,0) = 0.01 . The colours red, green and black correspond to the type-1 , type-2 and total population 
size respectively. D: A break of slope was observed in the population range. Initial propagules are »i (0,0) = 0.99,«2(0,0) = 0.01. (solid line) and 
«i(0,0) = 0.999,W2(0,0) = 0.001 (dashed line). 
doi:10.1371/journal.pone.0103409.g001 



standard deviation a wUl be referred to as the propagule diversity). 
Note that the parameter e'' gives the median dispersal ability of the 
initial propagules. Second, we group the individuals into n types, 
with the type i individuals having the dispersal abilities between 
(/— and id, witli 3 = {d^ — cf)/!!, and and cP' denote the 
maximal and minimal dispersal abilities, respectively. For simplic- 
ity, all individuals of type i were assigned the same dispersal ability 
df = ii-l)6 + S/2. 

The model was solved using the Fast Fourier Transform 
algorithm implemented in the SciPy library of Python [55]. To 
simplify the illustration, we call a type i individual a slower 
dispersal if df <d\ + {djj—d\)/'i, an intermediate disperser if 
dl+{dl-d^)/2<d}<dl + l(dl-d1)/l>, a fast disperser if 
dj >df + 2(d^—df)/3. The shape of the expected dispersal kernel 
k and the mean dispersal ability of the individuals at location x and 

n 

time t, ^ dfpi(x,t) were investigated. The instantaneous and 

;=1 



average rates of spread at time t were calculated as 
c,(t) = x*(t+l)-x*(t) and c(t) = (x''(t)-x*(0))/t, respectively, 
where x* (t) is the location of the range front defined for a certain 
threshold of detection u* as x*(t) = max{x>0,u{x,t)>u*}. Recall 
that we are interested in the spread of the total population u{x,t). 

Results 

Propagule with two dispersal abilities 

When the initial propagule contains two types of individuals, the 
number of the fast disperser type 2 individuals remained low at the 
initial phase, while the total population consisted mainly of type 1 
individuals, similar to the composition in the initial propagules 
(Fig. lA). The type 2 individuals gradually reached the expanding 
front through spatial sorting (Fig. IB) and then increased in 
numbers, with the type 1 individuals compressed back to the range 
core (comparing Fig. IB with Fig. IC). A breaking point of slope 
was observed during the range expansion (Fig. ID), indicating a 
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Figure 2. Asymptotic rate of spread of a population witli two dispersal abilities. Theoretical (black lines) and computed asymptotic rate of 
spread when two dispersal abilities are present in the population. Initial propagules are wi (0,0) = 0.99,M2(0,0) = 0.01 (blue triangle) and 
Mi(0,0) = 0.999,i/2(0,0) = 0.001 (magenta triangle). Other parameter values are d^ = l (1) r = 0.75 (2) r=l and (3) r = 2. A: Using Gaussian dispersal 
kernels. B: Using Laplace dispersal kernels. 
doi:1 0.1 371/journal.pone.01 03409.g002 



type II bi-phase expansion. The time lag to the breaking point 
decreased when more type 2 individuals were in the initial 
propagule but the asymptotic rate of spread was not afTected by 
the propagule composition (Fig. ID). 

The asymptotic rate of spread for type II bi-phase expansion 
(i.e. the rate of spread at the second phase of expansion) was 
heuristically derived as follows. We note that the solution of the 
model (Eq.(l)) becomes a travelling wave of the Uj{x,t) = e^'"''--^^^''^ 
which spreads at the following rate: 

c, = ^log{RM{di,Ai)) = mm\\og{RM(diJ.)), (3) 

where M(di,A) = [A:,(.s)e^'"'c& is the moment generating function 
of ki, and the interval / is (0,oo) for Gaussian kernels and is 
{f),y/l/dj) for Laplace kernels [7,56]. Recall that we are interested 
in the spread of the total population 
M = ui+M2 = e"''i*''"'i'' + e"^-2(-^-'-2'>. Because l2<h and 
X\C\=X2C2 for Gaussian and Laplace kernels, we have 
u{x,t)cce^^^^^^'^'^ as x approaches infinity, meaning that the 
spreading rate of the total population (Eq.(2)) is given by c = C2. An 
approximation has been derived by Lutscher [57]: 



(4) 



where is the kurtosis of the dispersal kernel ^2 ■ In particular, for 
Gaussian kernels we have: 



cK,^2\n(R)dl. 
For Laplace dispersal kernels we have: 



(5) 



2\n{R)dl 1 



(6) 



The derived rate of spread c fits well with the asymptotic rate of 
spread obtained from the numerical simulations with a Gaussian 
dispersal kernel and also sets a close upper bound for the rate of 
spread with a Laplace dispersal kernel (Fig. 2). 

Propagules with multiple dispersal abilities 

When a number of n dispersal abilities were present in the initial 
propagule, the process of spatial sorting gradually pushed fast 
dispersers to the advancing range front while compressing slow 
dispersers to the range core (Fig. 3). Spatial sorting was also 
detected by calculating the mean dispersal ability which kept 
increasing while expanding (Fig. 4A-B). The expected dispersal 
kernel has a fatter tail at the range front than at the range core 
(Fig. 4C-D); that is, individuals at the front are more likely to 
travel longer distances than individuals from the range core. 
Importandy, the expected dispersal kernel was found closely 
related to the dispersal kernel corresponding to the mean dispersal 
ability. Consequendy, for the expected Gaussian dispersal kernel 
can be approximated as. 



1 



■l2 



(x-yy 



(8) 



ilndHy,t) 

and the expected Laplace dispersal kernel can be approximated as. 
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Figure 3. A propagating population with many dispersal abilities. The parameter values are Uo = 100,K = 500,/i = 0.5,a = 0.25,c = 0.75. The 
colours red, green and blue correspond respectively to slow, intermediate and fast dispersers whereas the black lines represent the total population 
size at different time. 
doi:1 0.1 371/journal.pone.01 03409.g003 



1 



-\x-y 



2dHy,t) 



(9) 



where (f'{y,t) is the mean dispersal ability at location 31 and time I. 

We observed two distinct phases during the range expansion. 
The instantaneous rate of spread increased with time in the first 
phase (e.g. in the first 40 generations of the simulated population 
in Fig. 5), followed by a constant rate of spread in the second 
phase. Numerical simulations further suggested that the instanta- 
neous rate of spread can be estimated by the mean dispersal ability 
of the population at the advancing range front (Fig. 5). Specifi- 
cally, the instantaneous rate of spread c (t) can be approximated as 
the rate of spread of a single species with the kernel 



c*(f) = dl\n{K)d^(x*(t)), 



(10) 



where (f'ix'it)) is the mean dispersal ability at the front x'ii) of 
the invasion at time t. Furthermore, following the same procedure 
as for the scenario with two dispersal abilities (section 3.1), we 
found that both the instantaneous and average rate of spread 
approached a same asymptotic rate of spread: 



1. 



c= mm-\n{RM{d„),>yr 

X A 



21n(i^H2 1 



12 



(11) 



where is the maximal dispersal ability in the population. 

As the dispersal abilities of the individuals were randomly drawn 
from the lognormal distribution, In N(^,a), the maximal dispersal 
ability d^ is a random number. Let Z)(y„ be the random variable of 
the maximal dispersal ability in the initial propagxiles of size Uq. 



PLOS ONE I www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | el 03409 



Mixed Propagules and the Acceleration of Range Expansion 




10 
Location 



20 




20 
Location 



40 




Location 



100 



Location 



Figure 4. Evolution of the mean dispersal ability and expected dispersal kernels. (A-B) Mean dispersal ability: During the invasion, the 
dispersal ability at the front is higher than the dispersal ability at the core of the invasion. The region in blue indicates habitats that have not been 
invaded at each generation. (C-D): The expected dispersal kernel at the core of the invasion (in black) is narrower than that at the front of the 
invasion (blue). In both cases, the expected dispersal kernel (full circles) can be approximated by the kernel associated to the mean dispersal ability 
(Eq (8) and Eq (9)) (Solid lines). Figures (A,C) and (B,D) were obtained using Gaussian and Laplace dispersal kernels respectively. 
doi:1 0.1 371 /journal.pone.01 03409.g004 



The cumulative distribution function of Du^ can be given as, 
F{d,Uo) = P{Du,<d) = {^{d)f'', where ^(d) is the cumulative 
distribution function of the lognormal distribution, 
<^{d) = {\+erf{{\n{d)-li)/^/2a'-))|2, where erf{-) stands for the 
Gaussian error function. By solving F(d, Uo) =1/2 with respect to 
d we can obtain the median of Dy^: 



= exp 



(12) 



For Gaussian dispersal kernels with propagules having multiple 
dispersal abilities, the asymptotic rate of spread is, thus, given by: 



21og(i?)exp + 



1 



(13) 



This formula for the median asymptotic rate of spread was tested 
by solving Eq.(2) numerically for different propagule size Uq and 
propagule diversity ff. For each pair of f/o and (J the median of the 
asymptotic rate of spread from 15 simulations was calculated and 
compared with Eq.(13). Evidently, as the propagule diversity a 
approaches zero (i.e. all individuals have an equal dispersal ability) 
the asymptotic rate of spread becomes independent of the 



propagule size, whereas increasing either propagule size or 
propagule diversity will result in a faster spreading rate (Fig. 6). 

Discussion 

The role of propagule pressure in the introduction and 
establishment of introduced species has a rich background in 
literature [21-23]. In this work, we went further and investigated 
the contribution of propagule pressure in the actual rate of spread 
and the shape of the range expansion pattern. Two properties of 
the propagule were incorporated in the models, namely the 
propagule size (i.e. number of individuals in the propagule) and the 
composition of the propagule, namely the distribution of the 
dispersal abilities of the individuals in the initial propagule. Such 
propagule features were incorporated in our model using 
integrodifference equations (IDE). 

We considered two compositional structures of the initial 
propagule. In the fu'st case, we examined the importance of a 
small number of individuals with stronger dispersal ability in the 
propagule. A linear range expansion was observed during the 
initial phase of the invasion, which was followed by another linear 
expansion with a higher rate of spread. In other words, a sudden 
increase in the rate of spread can result from a mixed composition 
of the propagule. Furthermore, while the duration and rate of 
spread of the initial slow phase depended on the frequency of fast 
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Figure 5. Rate of spread of a mixed population. Parameter values 
are (/q = 100,_S: = 500,/( = 0.5,(7 = 0.25,)- = 0.75. Diamonds: Instantaneous 
rate of spread. Circles: Average rate of spread. Solid black line: Predicted 
asymptotic rate of spread. Solid red line: predicted instantaneous rate 
of spread (Eq (10)). The rate of spread predicted by a single population 
model is shown by the dashed line. 
doi:1 0.1 371 /journal.pone.01 03409.g005 

dispersers in the initial propagule, tlie long-term rate of spread of 
the population was determined by their dispersal abihty, however 
rare they were in the initial propagule (Eq (7)). This phenomena is 
expected when the fast dispersers do not go extinct, for example as 
result of demographic stochasticity which are important in 
populations at low density [58,59]. In the second case, we 
investigated the case where the dispersal abilities of the initial 
propagule are log-normally distributed. Dispersal abilities were 
subject to spatial sorting. More clearly, the expanding front was 
inhabited by individuals with higher dispersal ability than the core 
population, as observed in different real invasions [41,42,44]. 
Furthermore, the frontal population was less mixed (consisted only 
of fast dispersers), compared to the population at the core of the 
invasion where all dispersal abilities were present, as predicted by 
competition models in which dominant species win and limit 
weaker individuals' invasion [12,59]. Finally, spatial sorting was 
more apparent for more mixed propagules. 

The spatial sorting of dispersal abilities was reflected in the 
mean dispersal rate of the population. At the beginning of the 
invasion, individuals with better dispersal abilities are low in 
number as their growth is limited by the individuals with weaker 
dispersal abilities. However, fast dispersers do not go extinct as 
they are as competitive as the slow dispersers. As individuals with 
better dispersal ability reach the front of the invasion, they can 
grow in number in the open space without any competition against 
the slow dispersers. The mean dispersal ability at the frontal 
population therefore increases. This process occurs at every 
generation during range expansion, and results into an increasing 
dispersal ability on the expanding edge. The increasing dispersal 
rate in turn yields an increasing rate of spread, that is, an 
acceleration of the range expansion. Numerical simulations 
suggested that the rate of spread between two generations can 
be approximated using the dispersal ability of the frontal 
population only. This result is consistent to the findings of Bouin 
et al. [60] . This result suggests that ( 1 ) empirical quantifications of 



the dispersal ability, such as the mean dispersal rate for all 
individuals are only accurate for a short period of time and may 
underestimate the long-term rate of spread of the population 
(dashed line in Fig. 4) and (2) empirical predictions based on 
dispersal abilities obtained from the core population can depre- 
ciate the real rate of spread. 

Unlike other works which took possible mutations of dispersal 
relevant traits into account [8,61-63] the rate of spread obtained 
from our model remains bounded as the dispersal abilities are 
bounded. After the initial acceleration, the population expands at 
a constant rate. A close formula for the asymptotic rate of spread 
was derived for the constant asymptotic rate of range expansion 
(Eq. (14) with Eq. (13)). First, it is worth recalling that the rate of 
spread was obtained with the assumption that the dispersal ability 
of the propagule is log-normally distributed. This assumption was 
used due to different evidences that species-abundance relation- 
ships follow a log-normal distribution [53,64,65]. The rate of 
spread, however, can be derived for dififerent propagule cUstribu- 
tions simply by using the corresponding cumulative probability 
function. For instance, for a normally distributed initial propagule, 
the cumulative probability function is given by 
^(d) = {\ +erf((d — fi)/\/2(T^))/2 and the rate of spread is given 
by 



2\og(R){ fi + V2(j^erf-^ { 2 



j_ 



where Uo and a denote the propagule size and compositional 
diversity, respectively. 

Second, the obtained rate of spread is similar to the 
approximation for the RD model (c = 2i/rD). Moreover, the RD 
result is obtained when a common dispersal rate is shared by all 
individuals, by letting a tend to 0. Finally, the expression of the 
rate of spread suggests that the rate of spread increases with the 
propagule pressure. This result is in hne with the speculation that 
increasing the propagule size can improve the species spread by 
providing better suited individuals for invasion [22,66]. Further- 
more, our results are consistent with the findings of Skalski and 
Gilliam [30] and Yamamura [67] who explored different models 
to elucidate the importance of different dispersal abilities in a 
population. 

Despite the theoretical progress made here in understanding the 
acceleration of range expansion in biological invasions, our model 
only captured one facet of the dynamic nature of dispersal 
strategies, through the spatial sorting of individuals with diverse 
dispersal abilities in the initial propagule. Other factors do exists, 
particularly the evolutionary dynamics of dispersal-relevant traits, 
which can also affect spreading dynamics. Simulation models have 
suggested that the eco-evolutionary dynamics of dispersal can 
either promote [37,38] or prohibit [36,38] spatial sorting, and the 
acceleration of spreading especially in environment with high 
temporal variability of habitat [36] . Rapid evolution can further 
affect the ecological dynamics of a population during colonization 
and spread [68-70]. Further models therefore need to take 
evolutionary detail into consideration to improve the prediction of 
range dynamics under local and global environmental changes. 

To conclude, our results suggest that the variety of dispersal 
abilities in the initial propagule plays an important role in shaping 
the range versus time pattern during a population's spread. A 
biphasic invasion, which consists of two linear range expansions 
with different rates, resulted from a propagule with two dispersal 
levels. When the initial propagule was more mbced, the invasion 
started at a slow rate and then accelerated until a maximal rate of 
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Figure 6. Propagule pressure and the median rate of spread. Theoretical (left) and computed (right) median of the asymptotic rate of spread 
for r^O.75, A: //^l and B: //^2.5. The computed median was obtained from 15 simulations for each set of parameter values. 
doi:10.1371/journal.pone.0103409.g006 



spread was attained. In addition, our results emphasize the 
importance of census time and locations when estimating the 
parameters of reaction-dispersal models as data collected from the 
core of the invasion may underestimate the actual rate of spread. 
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